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Abstract 

The spatial development of an initially linear vorticity-mode instability on a compress- 
ible flat- plate boundary layer is considered. The analysis is done in the framework of the 
hypersonic limit where the free-stream Mach number M -¥ oo. Nonlinearity is shown to 
become important locally, in a thin critical layer, when <r, the deviation of the phase speed 
from unity, becomes o ^ Af — ^ ^ and the magnitude of the pressure fluctuations becomes 
O (cr iM 2 j . The unsteady flow outside the critical layer takes the form of a linear in- 
stability wave but with its amplitude completely determined by the nonlinear flow within 
the critical layer. The coupled set of equations which govern the critical-layer dynamics 
reflect a balance between spatial-evolution, (linear and nonlinear) convection and non- 
linear vorticity-generation terms. The numerical solution to these equations shows that 
nonlinear effects produce a dramatic reduction in the instability-wave amplitude. 

1. Introduction 

The theory of compressible boundary-layer instabilities has received renewed attention 
in recent years primarily due to the current interest in hypersonic flight. The foundations 
of the linear theory were provided by Lees & Lin (1946) who showed that compressible 
boundary layers are often inviscidly unstable due to the presence of a generalized inflection 
point. Indeed, at sufficiently large Mach numbers, it is believed that laminar boundary-layer 
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instabilities are predominantly inviscid. 

A comprehensive account of the linear inviscid stability theory is given by Mack (1984,1987). 
There it is shown that, in addition to the single inflectional mode present at all non-zero Mach 
numbers, the Rayleigh problem admits a countably infinite set of solutions whenever there 
is a region of supersonic flow relative to the phase speed of the disturbance. An example of 
this multiplicity of solutions is shown in figure 1 where the spatial growth rate — is plotted 
against the frequency S for an insulated flat plate at a free-stream Mach number M = 10. 

In this figure, Mack’s (1987) method of identifying families of eigensolutions associated with 
the neutral sonic modes is used. The first family, denoted E u , starts from the upstream 
running Mach wave at S = 0 where c = 1 — 1/M while the other families, denoted E^ n , start 
from the successive downstream running Mach waves at S = Sdn where c = 1 4- 1/M. Only 
solutions for c r < 1 are presented in the figure so the Ed n curves begin at the frequencies 
Sin corresponding to the non-inflectional neutral points with c — 1. The E^n curves are 
shown only for a small range of 5, were this range extended, the damping rate along each 
curve would, in general, reach a maximum and then decrease toward zero. The other neutral 
frequencies shown in figure 1 are denoted by S sn and correspond to the generalized inflection 
point. In this paper, the convention of referring to the E u family as the vorticity mode and 
the E dn families as the acoustic modes (with the first acoustic mode corresponding to E<n, 
the second to E^ and so on) will be followed. 

Single-frequency excitation of compressible boundary layers gives rise to spatially grow- 
ing instability waves. For excitation amplitudes and Reynolds numbers sufficiently small 
and large respectively, the initial evolution of such disturbances can be described by weakly 
nonparallel linear stability theory. According to this theory, the instability wave continues 
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to grow as it propagates downstream and, owing to the growth of the mean boundary layer, 
the local Strouhal number increases. The latter result follows from the definition of the local 
Strouhal number as SF* /Uq 0 where S is a measure of the local boundary-layer thickness, F 
is the constant forcing angular frequency and Uqo is the constant free-stream velocity. For 
simplicity, the mean flows considered here are assumed to be self-similar so the local stability 
properties can be found from eigenvalue diagrams such as figure 1 where 5 is now interpreted 
as the local Strouhal number. 

As the local Strouhal number increases towards its neutral value, nonlinearity will come 
into play locally, in a thin critical layer, to balance the resulting singularity in the Rayleigh 
equation provided the Reynolds number is sufficiently large (see Goldstein & Leib 1988). 
This nonlinear interaction can be described by an analysis completely parallel to that given, 
for compressible shear layers, by Goldstein & Leib (1989) and Leib (1991). For the shear 
layers, the near-neutral region is the first region in which nonlinearity enters the development 
of the unsteady flow provided the instability- wave amplitude remains small. However, for 
the boundary layers of interest here, nonlinearity can become important away from the near- 
neutral region while disturbance is still small due to the interplay between the vorticitv 
and acoustic modes. This was first noted in the large Mach number analysis of the Rayleigh 
problem done by Blackaby, Cowley & Hall (1992, hereinafter referred to as BCH). They show 
that near each non-inflectional neutral frequency Si„ the vorticity mode has a well-defined 
critical layer and consequently, may be preferentially affected by nonlinearity. Referring 
to figure 1, this means that if a vorticity mode, with a dimensionless frequency 0.1 say, is 
excited, nonlinearity can become important as the local Strouhal number approaches either 
5x2 or S 13 (far from the neutral Strouhal number S s4 ) while the instability-wave amplitude 
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is still small. It is this possibility that will be investigated in the present work. 

As in the compressible-shear-flow analysis of Goldstein h Leib (1989) and Leib (1991) 
and the incompressible boundary and shear layer analyses of Goldstein Durbin & Leib (1987), 
Goldstein &; Leib (1988) and Goldstein fe Hultgren (1988), nonlinearity comes into play 
locally, in a thin critical layer, to balance a singularity in the Rayleigh equation only now the 
singularity arises as the local Strouhal number approaches a non-inflectional neutral value. 
The unsteady flow outside the critical layer remains essentially linear but with a slowly 
varying amplitude that is completely determined by the nonlinear dynamics in the critical 
layer. The appropriate scaling between the local instability-wave amplitude and frequency 
for which this nonlinear interaction takes place is determined by requiring that the nonlinear 
problem match onto the linear small-growth-rate solution far upstream. This ensures that 
the nonlinear solution represents the natural downstream continuation of the upstream linear 
solution. 

To obtain a description of this nonlinear interaction from a first principles analysis, it is 
necessary to characterize, in some rational way, the relative sizes of the linear growth rate and 
wavenumber. This characterization is most easily arrived at by considering the problem in 
the hypersonic limit where the free-stream Mach number M — ► oo. The relevant asymptotic 
solutions to the Rayleigh problem, for a fluid obeying Chapman’s viscosity law, are then 
given by Cowley & Hall (1990) and Smith & Brown (1990). The former reference contains an 
analysis of the small-wavenumber acoustic modes while the latter gives order-one wavenumber 
solutions for the vorticity mode as well as the acoustic modes. A detailed analysis of both 
mode types is given by BCH for a fluid satisfying the more accurate Sutherland viscosity law 
and Grubin & Trigub (1992a, b) present solutions for a fluid with a Prandtl number less than 
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unity and a power law viscosity-temperature relation. 

The nonlinear interaction considered here is inherently a large Mach number phenomenon 
because it requires that the first few S ln lie within the frequency range of the unstable vorticity 
mode (c.f. figure 3). But, from a practical stand point, the Mach number should not be too 
large since as it increases the frequency of the most rapidly growing vorticity-mode solution 
(and consequently, the excitation frequency of most interest) moves beyond the range in 
which the interplay with the merged acoustic modes produces a well-defined critical layer 
(see BCH). For an insulated fiat plate, these considerations restrict the Mach number to 
values between about 8 and 20 which, nevertheless, is a range of technological interest. BCH 
indicate that their asymptotic solution for the vorticity mode is not in overall agreement 
with the exact solution until M > 20. It is felt that the nonlinear evolution equations to be 
derived here will remain valid in the more moderate range, 8 < M < 20, for the following 
reasons. The present analysis involves only the linear solutions for the vorticity mode near 
the non-inflectional neutral frequencies Si„ and these solutions are in, at least, qualitative 
agreement with the exact solutions at lower Mach numbers. Also, the nonlinear interaction 
is local in nature and, to a large extent, its occurrence depends only on the requirement that 
the linear growth rate be small compared to the wavenumber which is certainly the case in 
the Mach number range of interest. 

Goldstein & Wundrow (1990) showed that the aforementioned nonlinear mechanism 
causes a breakdown in the linear solution of Cowley & Hall (1990) when the amplitude of 
the pressure fluctuations in the main boundary layer becomes O (m - 4 In -1 M 2 ). Although 
presented in the context of acoustic modes, these analyses also apply to the vorticity mode 
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in the vicinity of the non-inflectional neutral frequencies 1 . The solutions to the nonlinear 
critical-layer problem derived in Goldstein & Wundrow (1990) depend on a single free para- 
meter, r, which is defined as the ratio of the logarithmic derivative of the mean temperature 
to the logarithmic derivative of the mean vorticity evaluated at the transverse position where 
the phase speed of the disturbance equals the mean flow speed. This parameter is between 
0 and 1 on the acoustic-mode branch of the linear solution and is greater than 1 on the 
vorticity-mode branch. The authors only presented numerical results for values of r less than 
1, however these results clearly show the significance of nonlinear vorticity generation due to 
compressible effects. In particular, for values of r near 1 they found that these effects come 
into play first to produce a dramatic increase in the instability-wave growth rate. 

The purpose of the present work is to determine how critical-layer nonlinearity affects 
the spatial development of an initially linear vorticity mode as it nears a non-inflectional 
neutral frequency. This objective could be achieved by computing solutions to the evolution 
equations derived in Goldstein & Wundrow (1990) for values of r greater than 1, however 
two changes will be made to the derivation of these equations which will greatly broaden 
their range of applicability. First, Sutherland’s law will be used to model the variation of 
the viscosity with temperature instead of Chapman’s law since, as pointed out by BCH, the 
former law is better suited for the large temperature variations encountered at hypersonic 
speeds; and second, the frequency range over which the nonlinear analysis applies will be 
enlarged to capture the maximum possible deviation from S i n for which the linear vorticity 

J The approximate linear dispersion relation obtained in Cowley Hall (1990) and Goldstein & Wundrow 
(1990) has two branches for which the growth rate is positive. One branch corresponds to the amplified 
acoustic modes and the other corresponds to those amplified vorticity-mode solutions lying near the non- 
inflectional neutral frequencies S i n . This issue is explored in some detail by BCH. 
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mode has a well-defined critical layer. 

In the present analysis, nonlinearity first becomes important when the amplitude of the 
pressure fluctuations in the main part of the boundary layer becomes O (a* where a is 
a measure of the deviation of the local phase speed from unity and must be o The 

coupled set of nonlinear equations that determine the critical-layer dynamics are a special 
case of those obtained by Goldstein &; Wundrow (1990) and must be solved numerically. The 
computations show that, rather than enhancing the growth rate as it did in Goldstein & Wun- 
drow (1990), nonlinearity now causes a dramatic reduction in the instability-wave amplitude 
- presumably because the nonlinear vorticity-generation (or Bjerknes) term appearing in the 
final critical-layer vorticity equation changes sign relative to the linear inhomogeneous term 
when solutions for the vorticity mode are considered. 

The overall plan of the paper is as follows. The problem is formulated in §2, where the 
nonlinear flow is shown to gradually evolve from the strictly linear large M solution. The 
unsteady flow outside the critical layer is a linear inviscid perturbation about the hypersonic 
(i.e. M > 1), non-interactive, boundary-layer flow analyzed in BCH and the solution for this 
flow is worked out in §3. In §4, a comparison between the asymptotic and exact solutions to 
the linear dispersion relation is made. The asymptotic solutions obtained in §§2 and 3 are 
then used to formulate the relevant critical-layer problem in §5 and the numerical solution to 
this problem is presented and discussed in §6. 

2. Formulation 

As in Goldstein & Wundrow (1990), concern here is with a nearly-inviscid, compressible, 
boundary-layer flow of a perfect gas over a flat plate. The flow parameters in the free stream 
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are used as reference quantities and are generally denoted by the subscript oo. The reference 
length Sq is taken to be some suitable boundary-layer thickness (e.g. momentum thickness) in 
the region where nonlinearity first becomes important. The steady flow is then characterized 
by the Mach number 

, ( 2 . 1 ) 

Goo 

and the Reynolds number 

P 1 t^oo^O 

1+C Voo ' 

where a <*> = y/yffifoo is the speed of sound in the free stream, C « 110.4/Too is the Sutherland 
constant for air temperatures measured in degrees Kelvin, v is the kinematic viscosity, 7 is the 
isentropic exponent of the gas and 3£ is the gas constant. When numerical values are needed, 
7 is taken to be 1.4 and C is taken to be 0.5 which corresponds to a typical free-stream 
temperature in the upper atmosphere. 

The Reynolds number is assumed to be large enough so that the unsteady motion is 
essentially inviscid and unaffected by mean boundary-layer growth over the length scale on 
which the nonlinear interaction takes place (see (5.13) for a more precise restriction). To 
the required level of approximation, the mean pressure is then constant and the local mean 
velocity and temperature are given by 

{uo,M = {/'(»?), 0}, (2-3) 



and 

0o = T{ 77), (2.4) 

respectively, where rj is the Dorodnitsyn-Howarth variable (Stewartson 1964) which is related 
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to the transverse coordinate y by 


y= r T(s)ds. (2.5) 

Jo 

The local mean velocity and temperature fields are, of course, dependent on the con- 
stitutive relations used to describe the fluid. The simplest fluid to consider from an analytic 
stand point is a perfect gas with a Prandtl number of unity and a viscosity-temperature rela- 
tion given by Chapman’s law since, in this case, the mean profiles can be expressed in terms 
of the Blasius function. However, the large Mach number behavior of both the steady and un- 
steady solutions can be significantly altered when more realistic models are used. The most 
dramatic modification results when the viscosity law is changed. If Sutherland’s viscosity law 
is used the steady flow takes on a double-layer structure in the limit as M — > oo (Freeman 
& Lam 1959) and BCH show that the asymptotic solutions to the Rayleigh problem proceed 
in inverse powers of M rather than \/InM as in the case of a Chapman-law fluid. A similar 
result is found by Grubin k Trigub ( 1992a, b) when the viscosity-temperature relation is given 
by a power law. The next most significant change occurs when Prandtl numbers other than 
unity are considered. This was investigated by Fu, Hall k Blackaby (1990) as well as Grubin 
k Trigub ( 1992a, b) and both show that the scaling of the mean flow remains unaltered but 
that the powers of M in the resulting asymptotic expansions depend on the Prandtl number. 
Finally, Fu, Hall k Blackaby (1990) considered the effects of gas dissociation and found that, 
despite modifying the quantitative results, these effects leave the qualitative results (i.e. the 
Mach number scalings and exponents) unaltered. 

As previously stated, the fluid in the present investigation is taken to be a perfect gas 
with a viscosity given by Sutherland’s law. It is now further assumed that the Prandtl number 
is unity. This yields a model that captures the essential large Mach number behavior needed 
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to formulate the nonlinear problem without resulting in an excessive amount of algebra. 
Accordingly, the function / in (2.3) is determined by 

//,, + 2 (f+c /# ') =0 (26) 

subject to the boundary conditions 

/( o) = /'( 0) = 0, /'(oo) = 1, (2.7) 

and the local mean temperature is expressed as 

T = T b + (1 - T b )f + ^M 2 (T 6 + /')( 1 - /') (2.8) 

with T b = 1 for an insulated wall. 

Freeman & Lam (1959) found that in the hypersonic limit the solution for / takes on a 
double-layer structure consisting of a wall layer where 

£ = M 27? (2.9) 

is order one and a temperature adjustment layer where 77 is order one. In the wall layer, / 
expands like 

/ = M-i/o(£) + ---. (2-10) 

Substituting (2.9) and (2.10) into (2.6)-(2.8) and equating like powers of M, one obtains at 
leading order 

/o/o+ 2 (^) =0 (2.11) 

with 

/o(0) = /o(0) = 0, /o = £ + £oo + ^ pH as ^ oo (2-12) 
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where 


T 0 = ^{T h + f 0 ) (1 -/o), 


b = 


72 


(2.13) 

(2.14) 


(7- i)(r 6 + i)’ 

and ^oo is a constant which must be determined numerically. The temperature in the wall 
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layer is 0{M 2 ) so it follows from (2.5) and (2.9) that this region has a thickness of 0{M 5) 
in the physical coordinate y. 

The algebraic approach to the free-stream conditions exhibited in the wall layer becomes 
exponential in the temperature adjustment layer where 


f = T} + M 2foo +M </2(j?) + 


(2.15) 


Substituting (2.15) into (2.6)-(2.8), g' 2 is found to satisfy 


Vd" + 2 


(l-ifr-lKIV+lW + C 9 ’/ 


= 0 


(2.16) 


subject to 

g'= — -^ + • • • as T)-> 0, g' 2 ->-0 as 77-^00. (2.17) 

T 

Since the temperature here is 0(1), the width of this region in terms of y is also 0(1). 
Therefore, in physical coordinates, the wall layer makes up the main part of the boundary 
layer with the temperature adjustment layer serving as a relatively thin region through which 


matching with the free stream is accomplished. 

The unsteady motion is assumed to start as a linear, inviscid, spatially growing wave up- 
stream in the flow. It is further supposed that the wave is monochromatic (at least initially) 
and two-dimensional since, in the Mach number range of interest, the most rapidly growing 
linear waves are planar (Mack 1984,1987). The nonlinear interaction to be considered here 
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occurs as the local Strouhal number approaches a non-inflectional neutral value Si n . Non- 
linearity then affects the evolution of the unsteady flow over a streamwise scale that is long 
compared to the wavelength of the disturbance but short compared to the scale on which 
the mean flow develops. For the purposes of the present analysis, S can then be taken as a 
constant and given by 

5 = -JJ— = Si n + a 2 MiSi as M oo (2.18) 

where Si n has an expansion in inverse powers of M, S\ is an order-one constant and ct 2 M 2 
1. The behaviors of the local phase speed c and complex wavenumber a depend on both the 
sign of 5i and the size of the small parameter <r. For flow over a flat plate at zero angle of 
attack, the local boundary-layer thickness and consequently, the local Strouhal number both 
increase monotonically with increasing downstream distance. The local Strouhal number 
then approaches 5i n from below so, in the streamwise region of interest, Si is negative. The 
size of a is chosen so that the linear wave has a well-defined critical layer. The analysis of 
BCH shows that this condition is satisfied if 

M ~ 2 < <r < M~% (2.19) 

in which case 

c — 1 - (7 c, H , (2.20) 

Rea = M _ 2 q 0 -+ , (2.21) 

and 

Ima = — cr~M*K r + - • •, (2.22) 

as M — ► oo, where ci, do and n r are all order-one real constants. It should be noted that 
the so-called acoustic-mode scaling, on which the analysis of Goldstein & Wundrow (1990) is 
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based, is recovered, for the present choice of viscosity law, when a — M ~ 2 . The analysis to 
follow, however, does not apply to this scaling and a separate treatment would be necessary 
to capture this smaller range of frequency deviations. 

The critical layer is centered about the transverse position where the phase speed of the 
disturbance equals the mean flow speed. From (2.3) and (2.20), this requires /' — 1 = 0{v) 
which, in view of (2.15), (2.17) and (2.19), is satisfied in an overlap region between the wall 
and temperature adjustment layers where 

z = <74 M2 77 (2.23) 

is order one. The flow will be nonlinear in this layer provided the local instability-wave 
amplitude is sufficiently large. 

The unsteady flow remains predominantly linear outside the critical layer so the relevant 
solutions for the streamwise and transverse components of velocity, the temperature and the 
pressure expand like 


U = f' + cRe'I'i AV* + e 2 u 2 + * M 4 u 3 H , 

(2.24) 

v = -M“2 tReia^i A f e iJ? + e 2 t> 2 + M 4 v 3 H , 

(2.25) 

9 = T + tRe©iA + e iJ? + e 2 6 2 + M 4 9 3 + 

(2.26) 

and 


py = 1 + tRelli AV* + e 2 n 2 + €<t^T M 4 jt 3 -i , 

(2.27) 

respectively, where 


X = M~*a [x - (1 - <rc) t] 

(2.28) 
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denotes the streamwise coordinate in a reference frame moving with the wave and e <C 1 is a 
characteristic amplitude of the instability wave in the region where nonlinearity first comes 
into play. This nonlinearity alters the solution over a streamwise scale inversely proportional 
to the local linear growth rate and is accounted for by the function ^(ari) where 

x ! = (2.29) 

The 0(e 2 ) terms in the above expansions account for higher-order corrections due to nonlin- 
earity outside the critical layer while the 0(ea ** M 4 ) terms are introduced to allow matching 
with the higher harmonics generated within the critical layer. 

The relationship between e, a and M is determined by requiring that the lowest-order 
nonlinear terms appearing in the critical-layer vorticity equation are of the same order as 
the linear convection and non-equilibrium terms. This choice ensures that the nonlinear 
problem will reduce to the appropriate linear solution far upstream. The function A^(ari) will 
ultimately be determined by the nonlinear dynamics in the critical layer but matching with 
the upstream linear solution requires 

->• c t e (Kr+i ' t<)l1 as x\ — ^ — oo (2.30) 

where at is a complex constant. 

The local instability- wave amplitude e is related to the amplitude of the imposed dis- 
turbance through the solution of the weakly nonparallel linear stability problem. For a given 
excitation amplitude, the e-a-M relation then determines the magnitude of the deviation of 
the local Strouhal number from 5i n in the streamwise region where nonlinearity is expected 
to become important. The success of the theory for any given excitation amplitude depends 
on how well (2.21) and (2.22) approximate the wavenumber in that region as determined 
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by the local compressible Rayleigh problem. A further constraint on the theory is given by 
(2.19) since this condition ultimately leads to a restriction on the size of e (see (5.15) below) 
and consequently, on the excitation amplitude. 

The functions III, $ 1 , and ©i of 77 , x u a and M are determined, to the required 

level of approximation, by 

(/' - cf 4- [ — ^—2^-1 - <* 2 r 2 fi - M 2 (// ~ c)2 ] n x = 0 , (2.31) 

U ’ dr, |_(/'- c ) 2 dr, \ [ T J 

4- (-P-) = ~T - ll n a , (2.32) 

dr, \f — c) L (f'-c y J 


^ = - M ' 2 jrr c Ul + 


©1 = (7 - 1) Till 4- Yjjt _ - ^ 1 ’ 


subject to the boundary conditions 


— — = 0 at 77 = 0, III -4 0 as 77 -4 00 , 


where the complex wavenumber and phase speed are now given by 


a = M 2 I a + <7 4 M 1 - 


1 d 


i At dxi / 1 


u 2 1 dA^ 
= 1 — crc - a * M — r-j — > 
iaA T dxi 


respectively, and the real constants a and c have expansions of the form 


Ot — OtQ + (T 2 M 2 Ot\ + 


c = ci H . 


(2.39) 
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It should be noted that the next correction to the mean flow in the wall layer was found by 
BCH to produce an O (m _ 5( 1+v ^ term in (2.38). Although the essential results for the 
linear flow outside the critical layer are unaltered by this term, (2.19) is stiffened slightly by 
requiring 

a > (2.40) 

so as to keep the linear analysis as straight-forward as possible. 

In the next section, approximations for n l5 4>x, \Pi and ©i valid outside the critical layer 
are derived. To facilitate in this derivation, a schematic of the asymptotic structure of the 
flow field in the Dorodnitsyn-Howarth variable is given in figure 2. 


3. Unsteady solution outside the critical layer 


First, the unsteady solution in the wall layer where £ = 0(1) is considered. In this 
region, IIi expands like 

III = P 0 + ct 2 M 2 Pi + • • • + i a^M 2 Pi + ■■■ (3.1) 


where, without loss of generality, Pq is taken to be real so that i Pi is the first imaginary term. 
Substituting this expansion, along with (2.9) and (2.10), into (2.8) and (2.31) and equating 
like powers of a and M the following equations are obtained 


CPq = 0, 


(3.2) 


and 


CP i 



o 2 0 T 2 



(3.3) 


CPi = 21m 




OS-1)/ 5 1 

To 



(3.4) 
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where 


£s(/i “ 1)J l[urT?l | ' 4r ' 11 ' 


(/o - ir 

To 


In the frequency range of interest, the eigenfunction is concentrated in the wall layer so, in 
addition to the usual wall boundary conditions 

i*(0) = P[(0) = J?(0) = 0, (3.6) 

the solutions to (3.2)-(3.4) must also satisfy 

P l ~2^B a C l , P, ~ -21m B » C „ ( 3 7 ) 

as ^ oo, where Bo is a nonzero constant which depends on the normalization chosen for 
III and the constants Ci and C t are found by the method of variation of parameters to be 

and 

c - = mI:( j^)’ (**+«**)«. (3 ' 9) 

A numerical solution to (3.2) subject to (3.6) and (3.7) yields a countably infinite set of 
eigenvalues cio each of which corresponds to a non-inflectional neutral frequency S\n- 

The solution given by (3.1) breaks down in the overlap region where z = 0(1). The 


appropriate expansion, in this region, is 


III = a 2 M 2 Eo + • — I - cm Pi + h i<T 2 M 2 Pi + • ■ • 


(3.10) 


where E 0 is a real constant and Pi and iP t contain the first non-constant real and imaginary 
terms respectively. It is easy to see from (2.9), (2.10), (2.12) and (2.23) that in this region 


/ = 1 _ CT - + 


(3.11) 
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Substituting (3.10) and (3.11), along with (2.23), into (2.8) and (2.31) one finds, after integ- 
rating twice with respect to 2 , 


- „ / b 2 2bci , \ „ 

P ' =D ‘ -7? + 3^ + ^ 2 +E *' 


(3.12) 


and 


«=- 35 f 


+ 


2 bci 

Yz J 


C,z 


-|-2Im 


' 1 d A'\ _ / b 

iYYY; Di 


) + Ei 


(3.13) 


where D\, E\, Di, Ei and 4> ± are at most functions of X\ (the latter of which can be different 
depending on whether z ^ (6/ci)< ). 

The eigenfunction must be rescaled again in the temperature adjustment layer. The 
solution here turns out to be 


IIi = ct 2 M 2 Eq + h (7 2 M2 [F\ij + G \ ) -I f \a * M 2 ( Fit) + <?,) + 


(3.14) 


where only the first terms containing non-constant real and imaginary parts are shown. 
Clearly, (3.14) does not satisfy the free-stream boundary condition given by (2.35) and it is, 
therefore, necessary to introduce an outer region where the variable 


C = M 27 ? 


(3.15) 


is order one. The solution in this region is given, to the required order of accuracy, by 


IIi ~ Q ex P (— oiy/l — o 2 M 2 c 2 


(3.16) 


where the constant Q has an expansion in powers of a and M. 
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The solutions in the various regions must, of course, be matched. The relevant results 


of this process are 


and 


Eq = 2-^-BoCi, 

(3.17) 

0 0 


Di = Bo , 

(3.18) 

D, = -Im</T, 

(3.19) 

- c\ 

1 “ 2Cj ’ 

(3.20) 

= "bal 4 Im (<f> + - 4 >~) , 

(3.21) 


where the last two equations form the leading-order approximation to the linear dispersion 
relation. The phase shift <t> + ~4>~ appearing in (3.21) is ultimately determined by the nonlinear 
dynamics inside the critical layer but matching with the upstream linear solution requires 


<f> + — <f> — > nr as Xy — > — oo. 


(3.22) 


In order to formulate the nonlinear problem in the critical layer, it is necessary to have 
the expansions for $ ls and ©i near the transverse position where the phase speed of the 
disturbance equals the mean flow speed. These are obtained by substituting (3.10)-(3.13) 
along with (2.23) into (2.32)-(2.34) and expanding about z = z c = ( b/ci )«. The resulting 
expressions are 


.,34ft, n J5o ii l / OoCi 1 d^^ Bo 

*■ - (324-^3- + — _ j ^ 


+ 


1 b » B o , --2 

*1 = z )-=5 ~ c M 


9 z| ' ' 05 


Jo + -r- (In | z - z c \ + 4^ - <t> )| cj ^ + 
ft J <*0 


(3.23) 

(3.24) 
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and 


©i + (3.25) 

z c <* 0 

where the real constant B has an expansion of the form 

B — Bq -h • • • + o (^7* M 2 ^ , (3.26) 

Jo is an order-one real constant and the dots indicate that higher-order terms in z — z c as well 
as in a and M have been omitted. It is interesting to note that the leading-order temperature 
fluctuation given by (3.25) is bounded in the critical layer. This result, which differs from the 
singular behavior obtained by Reshotko (I960), is due to the long wavelength approximation 
being employed here and is the primary reason the critical-layer dynamics turn out to be 
strongly nonlinear rather than weakly nonlinear as in the compressible-shear-flow analysis of 
Goldstein & Leib (1989) and Leib (1991). 

4, Linear dispersion relation 

Before proceeding to a formulation of the nonlinear problem, some results obtained from 
the linear analysis are discussed. 

The local Strouhal number of the instability wave can be expressed as 

5 = ac = M~i do + <J 2 M^6t\ H , (4.1) 

so it follows from (2.18) that 

Sin = M~td 0 + o ( <t 2 M 2 ) , (4.2) 

and §1 = di. Values of (4.2) corresponding to the first six do determined by solving (3.2) 
subject to (3.6) and (3.7) are shown as a function of M in figure 3 for an insulated plate. 
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Also shown are the exact values of S\ n determined by solving (2.31) subject to (2.35) with 
c = 1. This figure demonstrates that, even at supersonic Mach numbers, the leading-order 
hypersonic solution provides a good approximation to Si n . 

In addition to the 5i„ curves, the exact inflectional neutral frequencies S sn are shown in 
figure 3. These values of 5 correspond to c = f'(r) s ) where T) s is determined by the generalized 
inflection point criterion, 

Y (j^) = 0 at r} = T)s - (4 - 3) 

As first observed by Mack (1984,1987) , each S sn curve has a range of M for which d S/d M > 0 
and the corresponding neutral mode is vortical in character. In general, the rising portions 
of the curves mark the high frequency neutral point of the vorticity mode while the declining 
portions mark that point for the acoustic modes (c.f. figure 1). Figure 3 suggests that 
the slope of the rising portions tends to zero as the Mach number increases in contrast to 
the behavior exhibited in figure 6(a) of Cowely k Hall (1990) which was computed for a 
Chapman-law fluid. Smith k Brown (1990) found that along the rising segments in the latter 
figure S is 0(V In M) asM^oo while, for the Sutherland-law fluid considered here, BCH 
show that 5 is 0(1). 

A direct comparison with Mack’s (1984,1987) work is difficult since he used a combin- 
ation of viscosity laws and took the Prandtl number to be 0.71. Still, the neutral curves 
in figure 3 as well as the growth rate curves in figure 1 are in qualitative agreement with 
his results. The assumption of a Prandtl number of unity made here appears to yield solu- 
tions which reasonably approximate those computed with a more accurate value. However, 
it should be noted that Grubin k Trigub (1992b) found significant changes to occur in the 
Rayleigh solutions when the Prandtl number is sufficiently small. 
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Having demonstrated an agreement between the exact and asymptotic properties of the 
neutral solutions, the behavior of the amplified waves is now considered. The local complex 
wavenumber of the linear disturbance is given by 

a = j^d 0 + <7 2 M 2 oti H ) + i (-a” M 2 k t H ) j (4.4) 

where it follows from (2.30) and (3.20)-(3.22) that 

K r = 2» • 3 5 7rA (-di) T (4.5) 

with 

A = . (4.6) 

b* 

Equations (2.14) and (3.8) show that b and C\ are real and positive, consequently, (4.5) 
requires negative values of di to yield real values of n r . When di is positive, the growth rate 
was found by BCH to be 0(crM~i) so it follows from (4.1) and (4.2) that a discontinuity 
in the expansion for Ima arises as S passes through 5i„. BCH went on to show how this 
discontinuity is smoothed out by the solution for the acoustic-mode scaling a — M~ 2 . This 
aspect of the problem need not be considered here however, since, as mentioned in §2, Si 
(= di) is strictly negative in the streamwise region of interest. 

The variation of K r with non-inflectional neutral frequency and wall temperature is de- 
termined by A. This quantity is shown as a function of Tb in figure 4 for the first two 
values of do- The curves show that, for fixed values of Tb, di and M, the growth rate of 
the vorticity mode near the first non-inflectional neutral frequency Sn is less than that near 
the second S 12 which is in agreement with the results shown in figure 1. The asymptotic 
solution predicts that this trend of increasing growth rates is continued for the subsequent 
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Sin, yet, in figure 1, the growth rate of the vorticity mode near S \3 is clearly less than that 
near Si 2 * This discrepancy is due to the fact that at M = 10 the third acoustic mode has 
only just merged with the vorticity mode (c.f. figure 3). At M - 20, figure 6.2 of BCH 
shows that the growth rate near S 13 becomes greater than that near 5i2 in accordance with 
the asymptotic solution. If the growth rate associated with a particular non-infiectional neut- 
ral frequency is considered, figure 4 shows a decrease with decreasing T This result may 
seem to contradict the conclusion that wall cooling has a destabilizing effect on the Mack s 
so-called higher modes (Mack 1984). However, figure 9.13 of Mack (1984) clearly shows that 
while the magnitudes of the peaks in the vorticity-mode growth rate increase with decreasing 
wall temperature, the growth rate in the valleys actually decreases. 

Despite this qualitative agreement, quantitative agreement between the exact and asymp- 
totic growth rates is not achieved until the Mach number becomes very large. Quantitative 
agreement is, however, not crucial in the present investigation because the final nonlinear 
problem turns out to be entirely independent of K r . In fact, the nonlinear dynamics depend 
only on the a-M scaling of the growth rate. To determine the accuracy of this scaling, the 
exact vorticity-mode growth rate is computed as a function of M for a fixed deviation from 
a non-inflectional neutral frequency, that is for 

5 = (1 - A )Sm ( 4 * 7 ) 

where 0 < A < 1 is a constant. The results of this computation are presented in figure 5 
for the first two values of S ln . The asymptote predicted by the hypersonic solution is also 
shown. This is determined by noting that (4.1), (4.2) and (4.7) imply 

di ~ —a~ 2 M~ 2 Xdo (4-8) 


23 



which, from (2.19) and (2.40), requires 


1+V?) < A< M -f 


(4.9) 


Substituting (4.8) into (4.5) and the result into (4.4) gives 

a, = Ima = O {m~* A "S') 


(4.10) 


and, since A is held fixed along the computed curves, 

M da, 9 

rss , 

a, d M 4 

Although the comparisons shown in figure 5 are fairly good, it is felt that improvements could 
be obtained by adding higher-order corrections to (4.10). However, such terms do not enter 
into the critical-layer analysis and so are not pursued here. 


5. The critical layer 


Equation (3.24) shows that the solution in the overlap region becomes singular where 
z = z c . The governing equations must therefore be rescaled to obtain a bounded solution in 
this region. The thickness of the small-growth-rate critical layer is of the order of the mean 
velocity times the growth rate divided by the mean velocity gradient times the real part of 
the wavenumber so it follows from (2.36) and (3.11) that the appropriately scaled transverse 
coordinate in this region is 

Z = c- 7 *M~ 2 (z- z c ). (5.1) 

Equations (2.24)-(2.27), (3.10), (3.11), (3. 23)- (3.25), (A1)-(A4), (A 13) and (A 15)- 


(A 17) along with (5.1) suggest that the flow in this region should expand like 


u = 1 - crc-h ff l * M 2 ^rZ - cr*M 4 ^-Z 2 - (T*e-^-^ReA^e' x 
z\ z% 9z c a l 
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+ bounded Z — independent terms + <r 2 M 2 tu\ H , 


(5.2) 


v = a « M 2 e Re 


U b„ 324 _ 2 1 d A* 

75 Z + + rr- 


\OioA t dxi 


9z c ag / <*o 


+ 


,36 


144 


e = aM 2 ^ -a*M 4 ^Z + a*M 2 e0 1 + ---, 


and 


o'y =1 — cr 2 M 2 €^i5oReA^e i ‘^ + • • * . 


a o 


This solution will match the ‘outer’ solution in the overlap region provided 

(A -> * ^ReA'e'*, 

9z? a§ 

0i ->• — ^ReAV*, 

as Z ^ ±oo, and 

1 ,a- /■-“(' «5l _ » ^F teyl VS') e -« dZdJf = _^l BoW +_,i-)At. 

jt 7o J- co V dZ 9*f a 2 0 J b a 0 

The functions u x and ^ of I, Z and are determined by the vorticity and 
equations which can be written as 


u> - 3 7 1 [ ( _ dp n . -2 dp \ do 

D “ = -Dr + ^ M '^ p K“ajf + ) az 


dp ( dO ji 2 dO 

+ M ir 

dX dx 


_ djp_ ( 
dZ V 


;)]■ 


and 


respectively, where 


- 7—10- 

= ^ Dp , 

7 P 


ii 0 <9 5 3 v d 

D = <y~M u^+a{u-l + ac)-^+G , 


(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 

(5.8) 
energy 

(5.9) 

(5.10) 

(5.11) 
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(5.12) 


u> 


, ,_3 ( _ dv ji o <9u \ _3 , ,_3 1 du 

sM ‘{ a ax + c,M aTj-'’ ,M 1 t»z ’ 


and it has been assumed that 


i? > <r * M 2 


(5.13) 


so, to the required level of approximation, terms accounting for viscous diffusion effects as 
well as those arising from corrections to the mean flow due to the growth of the boundary 
layer can be neglected. 

The relationship between the amplitude scale €, the phase speed deviation a and the 
Mach number M is determined from (5.9)-(5.11) by requiring that the lowest-order velocity 
jump produced by the nonlinear terms is of the same order as the velocity jump due to linear 
effects, i.e. 0(a 2 M 2 €). A little experimentation shows that nonlinear effects will influence u\ 
if 


5 o 

£ = CT2M 2 


(5.14) 


which, when combined with (2.19) and (2.40), requires 

M -i(9+5x/7) < € < M -f (5.15) 


The equations governing u\ and Q\ are then 

® ($ - P) = ^ M Jz (*' - f • *) 


and 


where 


*(*-?*)- 0 

- d (A b„ b B 0 „ „ + :x\ d z* B 0n \(4b„ 

V = — (- o 0 ( t Z - — — r^ReA'e x ) — — + - 7 :^— Re ( —=Z 

dxi 0 \zf 9z c / dX 36 ao LV-zf 
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b z 3 a ° c i iaoAt dii 9 z c do 




(5.16) 


(5.17) 


(5.18) 
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is the limiting form of the total derivative in the critical layer. 

The presence of in the coefficients of both the % and d z terms in (5.18) as well as in 
the right-hand-side of (5.16) shows that nonlinearity enters the critical layer dynamics through 
both streamwise and transverse convection terms as well as through vorticity production due 
to pressure-density variations. The nonlinear streamwise convection terms turn out to be 
passive in that they do not affect the instability-wave growth rate (although they do strongly 
affect the flow field in the critical layer). These terms can be eliminated from (5.18) by simply 
introducing a strained transverse coordinate, 


Z = Z- ll^ReAV* 
36 c*o 

(5.19) 

so that V becomes 



(5.20) 

The parameters characterizing the ‘outer’ linear flow can be removed from the critical- 

layer equations by introducing the following normalized variables, 


o4 r 10 

A = -^%B 0 A'S x °- KiX '\ 
3 8 a 3 0 

(5.21) 


(5.22) 


(5.23) 

3 5 b _ 2 

X = — 0 : 0*1 *0, 

z z c 

(5.24) 

X = X — Xo + AC,' X 1 , 

(5.25) 

and 


y_2 3 Z< 6 {* + d*) 
~ 3 s So «»./' 

(5.26) 
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Equations (5.16), (5.17), (5.8) and (3.21) then become 



VQ = Re (i Ae' x ) ( 4 ^ - 3 ) , 

(5.27) 


VH = -Re (iAe 1 *) , 

(5.28) 


1 r2i r r+oo 

- / / ne~ iX dYdX = -3(<A + - <t>~)A, 

Tf 0 «/ — 00 

(5.29) 

and 

- r) = -to (^ ) . 

(5.30) 

respectively, where 

v=-T^W-^ X )w 

(5.31) 

If the, as yet, unspecified real constants Xq and xq are chosen to be 



X 0 = -arg (-.Boa 1 ) , 

(5.32) 

and 

1 /o4 -10 \ 

Io = -; to 

(5.33) 

(2.30) becomes 

A — * as x — > — oo. 

(5.34) 

From this condition along with (3.22) and (5.27)-(5.29), one can show that the solution for 

A remains real for all values of x. Equations (5.29) and (5.30) can then be combined to give 


1 r 2 * r+<*> . v dA 

- / Qe-^dydX = — i3— . 

7T JO J — 00 dx 

(5.35) 


It should be noted that equations (5.27) and (5.28) are the same as equations (6.29) and 
(6.30) of Goldstein &: Wundrow (1990) but with r = 4. The difference in the velocity jump 
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condition given by (5.35) and that given by equation (6.31) of Goldstein & Wundrow (1990) 
with r = 4 is due to a difference in normalizations which can be eliminated by replacing their 
T with (1 - r)T in equations (6.22)-(6.27) of the latter analysis. 

Equations (5.27), (5.28) and (5.35) must be solved numerically subject to matching 
with both the ‘outer’ flow in the overlap region and the linear flow upstream. The former 
condition requires that Q and H be strictly periodic in X and satisfy the homogeneous 
boundary conditions 

Q, H ->• 0 as Y -+ ±oo, (5.36) 

while the latter, in addition to (5.34), requires 

Q — 3 Re ( — (5.37) 

VT - ITT ) 

and 

H -+ -Re i ( 5 - 38 ) 

as x -4 -oo. The technique used to solve the final critical-layer problem parallels that given 
in Goldstein & Wundrow (1990) and the reader is referred there for a detailed discussion of 
the method. 

6. Numerical results and discussion 

Goldstein k Wundrow (1990) considered the spatial development of an initially linear 
inviscid instability wave on a hypersonic boundary layer. Their analysis was done for a fluid 
satisfying Chapman’s viscosity law using the acoustic-mode scaling which, in the present 
analysis, corresponds to o = M ~ The evolution equations they derived depend on a single 
free parameter, r, which is defined as the ratio of the logarithmic derivative of the mean 
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temperature to the logarithmic derivative of the mean vorticity evaluated at the critical level. 
For the amplified acoustic modes, r varies between 0 at the non-inflectional neutral frequency 
Sin and 1 at the inflectional neutral frequency S sn . For the amplified vorticity mode, r must 
be greater than some minimum value which is determined by the asymptotic solution to the 
Rayleigh problem and is itself greater than 1. This condition reflects the fact that the phase 
speeds of the vorticity-mode solutions to which the nonlinear theory applies are less than 
those of the amplified acoustic modes and so the critical layer is located closer to the wall 
where the mean temperature variations out weigh the mean vorticity variations. 

Goldstein k Wundrow (1990) only presented numerical results for the case when the 
disturbance is an acoustic mode, i.e. for r < 1. These results show that nonlinear vorticity 
generation in the critical layer can cause a super-exponential growth of the instability wave but 
that this subsides once the wave amplitude becomes sufficiently large. Transverse convection 
effects then come into play to produce a nonlinear roll-up of the flow within the critical layer 
and this, in turn, drives the growth rate toward zero. The final asymptotic state of the 
critical layer, for r < 1, is left undetermined but the authors indicate that viscous effects will 
eventually become important due to the continually decreasing scales generated by the roll- 
up. The incompressible-shear-flow analysis of Goldstein k Hultgren (1988) along with the 
investigation of the upper-branch stability of compressible boundary layers due to Gajjar k 
Cole (1989) suggest that viscous effects will ultimately produce a quasi-equilibrium critical- 
layer structure and a slow algebraic growth of the instability wave. 

Of primary interest, in the present investigation, is how the numerical results given in 
Goldstein k Wundrow (1990) are altered when r = 4, i.e. for the case of more practical 
importance when the disturbance is a vorticity mode. In figure 6, the scaled instability- 


30 



wave amplitude A determined by (5.27), (5.28) and (5.35) is shown as a function of the long 
streamwise variable x. The corresponding instability- wave growth rate is presented in figure 
7. These curves show that the wave grows linearly until its amplitude becomes large enough 
for nonlinear effects to become important. The growth rate then decreases and soon becomes 
negative. This results in a rapid reduction in the wave amplitude which is in contrast to the 
initial nonlinear amplification observed by Goldstein & Wundrow (1990). Both behaviors are 
the consequence of vorticity generation due to compressible effects. The difference presum- 
ably results because the nonlinear vorticity-generation (or Bjerknes) term appearing in the 
final critical-layer vorticity equation changes sign relative to the linear inhomogeneous term 
when r is greater than 1 (see equation (6.29) of Goldstein & Wundrow (1990)). The results 
presented in figures 5 and 6 correspond more closely to the special case found in the analysis 
of the weakly-nonlinear critical layer associated with a three-dimensional disturbance on a 
compressible shear layer done by Goldstein & Leib (1989) and Leib (1991). For this special 
case, the scaled amplitude function is purely real and figures 4 and 7 of the former analysis 
and figure 1 of the latter show that nonlinearity produces a rapid reduction in both the growth 
rate and wave amplitude. 

Although the nonlinear terms in the critical-layer problem become less important as the 
instability-wave amplitude continues to decrease, the solution does not return to its upstream 
linear state. This is because the effects of the nonlinear redistribution of vorticity, which takes 
place near the peak in the amplitude, persist downstream in such a way that the growth rate 
remains negative. This redistribution is demonstrated in figure 8 where constant-vorticity 
lines in the ( X , Y ) plane are shown at various values of x. As in Goldstein & Wundrow (1990) , 
localized regions of high vorticity are produced by the nonlinear coupling of the vorticity and 
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energy equations (c.f. figure 8(b)) while the more usual roll-up of the contours is the result of 
transverse convection effects (c.f. figure 8(c)). Superposed on this behavior is a straining of 
the contours due to the nonlinear streamwise convection term which has been re-introduced 


by showing the vorticity as a function of the more physical transverse coordinate, 


'-Sib 


fe*.) = Y -- 
4b do ) 8 


ReAe' x . 


(6.1) 


The straining is most pronounced in figure 8(c) where A is the largest. The important point 
to note is that, once the redistribution takes place, the vorticity field remains rolled-up even 
after the instability-wave amplitude has become quite small (c.f. figure 8(d)). 

The results in figures 6 and 7 suggest that the quasi-equilibrium critical-layer structure 
proposed for the downstream behavior of the acoustic modes will not apply to the vorticity 
mode considered here. In fact, the numerical solution to the critical-layer problem predicts 
that the scaled amplitude of the vorticity mode goes to zero at a finite downstream distance. 
Before this point is reached, however, the asymptotic formulation outside the critical layer 
breaks down. To show this, the amplitudes of the higher harmonics generated within the 
critical layer must be determined. The equation governing these amplitudes is obtained by 
equating the expression for the velocity jump across the critical layer given by the 'outer’ 
solution in the overlap region, (B 15), with that given by the 'inner’ critical-layer solution and 
is expressed as 


o4 r 17 i />27 t r+co 

ai m) = = -- / / Qe' mX dYdX (6.2) 

3 3 12 motlc\ 3 n Jo J - oo 

where m = 2,3, The amplitudes of the first few harmonics are shown in figure 9. It is 

clear that, as the amplitude of the fundamental goes to zero, the amplitudes of the higher 

harmonics do not. Consequently, the expansions (2.24)-(2.27) eventually become disordered. 
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A solution for the next stage of evolution is not worked out here because, once the amplitude 
of the fundamental becomes sufficiently small, other disturbances not accounted for by the 
present theory but always present in practice would then, most likely, be dominant. 

Finally, viscous effects could be incorporated into the critical-layer analysis by loosening 
the restriction in (5.13). Some care must be taken in doing this, however, since it is possible 
that an interaction between the mean flow and the leading-edge shock of the sort considered 
by Bush &; Cross (1967) could take place for certain values of a. In any event, it is felt 
that, although viscosity will tend to reduce the gradients in the vorticity and temperature 
fields, it will not significantly alter the behavior of the instability-wave amplitude because the 
nonlinear interaction that drives the amplitude to zero occurs over such a short distance (on 
the x scale) that viscous effects would not have time to assert themselves. 

Appendix A. Nonlinear terms in the overlap region 

The lowest-order approximations (in a, M and z—z c ) of the O (e 2 ) terms in the expansions 
(2.24)- (2.27) are determined in this appendix. Since these terms are generated by nonlinear 
interactions that are quadratic in the fundamental mode, they can be expressed as 


u 2 + n\U\ = ^ 2 °^ + Re'Pj^e' 2 *) 

(Al) 

v 2 + ffiVi = Rei^ 2 2 ^ e * 2 ^) > 

(A 2) 

9 2 + *!$! = G ( 2 0) 4- Re0^ 2) e i2 *, 

(A3) 

7r 2 = n' 0) + ReII< 2) e i2 *, 

(A 4) 


where «i, Vi, 0i and ttj denote the coefficients of the 0(e) terms in (2.24)-(2.27) respectively. 
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The functions II^, an ^ °f £i, a and M are determined to the required 


level of approximation by 


where 


L 2 n<’> = M 2 { a’rjf ’ - (/' - c) 2 ± [^j] } , 


$( 2 ) = I M~ 2 — - 

2 4a 2 (/' — c) l dr] 


Mi' Ml) 


*?* = jrr- c (y* 2> - M ~ 2 Tn i > ~ j*?’) . 
e| 2) = jrr- c [^*i 2) + (t - 1) (/' - c) rn< 2 > - iff” 1 ] , 


d I i d 


Xf'-c ) 2 


L » - v - <>' * - (Mr 2 [i - m 2 ^] , 

F 2 (2) = [2*? ($x*i) + A t2 , 

= ^©i ($1 ©1) - (7 - 1 ) (/' - c) IIx (©1 + Till) 

-( 7 - 1)T [2*iIIi (#xn r )] } At 2 . 


(A 10) 


(All) 


(A 12) 


Similar equations can be obtained for II 2° \ $2°^, ^2°^ an< t ©2^ however, these terms due not 
contribute to the leading-order behavior of 7 T 2 , 1)2, u 2 or 62 near the critical layer and so can 
be ignored. 

For the purposes of the present analysis, it is only necessary to consider the solutions in 
the overlap region where 2 = 0 ( 1 ). In this region, 11 ^ expands like 


i 4 2) = 4 2) + ---+*p{ 2) + 


(A 13) 
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where E ^ is at most a function of x x and Pf 2} is the first ^-dependent term. Substituting 
(A 13) along with (2.23), (3.10), (3.11) and (3.23)-(3.25) into (A 5) one finds, after integrating 


once with respect to z, 

dP | (2) 
dz 


- z c ) [z 2 c + 0(z - Z e )] (^ At ) 


to lowest order in a and M as 2 -> z c . It now follows from (A6)-(A8) that 

[z- 1 +0(z — z c )] (^^ f ) - 

*< 2 ) = o-h [J x + 0(z - z c )] (§A*V , 

\ a 0 ' 


and 


e (2> = [4 + 0(, - *,)] (fp) , 

as z — > z c , where J\ is an order-one constant. 


(A 14) 


(A 15) 
(A 16) 


(A 17) 


Appendix B. Higher harmonics in the overlap region 


In this appendix, the leading-order approximation to the velocity jump across the critical 
layer corresponding to the O {o 4 Af^e) harmonics in (2.24)-(2.27) is determined. Since these 
harmonics are generated by nonlinear interactions in the critical layer, they take the form 



*3 = Y, 
m= 2 

(Bl) 


00 

U3 = —M~i ^ Reimo^^d^e""*, 

m— 2 

(B 2) 


00 

0, = £ Ree<” ) 4 " ) e'“-', 

m =2 

(B 3) 

and 

» 3 =f;R*nWa<”'e“, 

771 — 2 

(B4) 
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where the scaled harmonic amplitudes a^^xi) are determined by the nonlinear flow within 
the critical layer and the functions ©3”^ and of 77, 21, a and M are 

determined by the following linear equations outside the critical layer, 


L m n| m) = o, 

(B5) 

(m) _ m~ 2 d ( TO ) 

3 m 2 a 2 (/' — c) drj 3 

(B 6) 

«4 m) = j±r c - M- ! rn‘"’) , 

(B 7) 


and 

0 3 m) = jr=~ c [f^ m) + (t - 1) (/' - c) rn< m) ] . (b 8 ) 

These functions satisfy the usual homogenous boundary conditions at the wall and decay 
exponentially as 77 — 00. However, they are, in general, discontinuous across the critical 
layer. 

It is only necessary to consider the solutions in the overlap region where 2 = 0(1). 
In order for the solutions in this region to match with those in the wall and temperature 
adjustment layers, II^ must expand like 

l4 m) = l + + + (B 9) 

for z < z c , and 

l4 m) = 1 + • • • + <7-*M~ 2 p[ m)+ + • • (B 10) 

for z > z c , where is the first ^-dependent term. Equations (B9) and (B 10) imply that 

the pressure fluctuations are constant across the critical layer to 0{<j~ M 4 e) (a fact that can 
be deduced from the transverse component of the momentum equation in the critical layer), 
but that the pressure gradient is discontinuous at 0(cr2 M 2 e). Substituting (B 10) along with 
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(2.23), (3.10), (3.11) and (3.23)-(3.25) into (B5) one finds, after integrating once with respect 


to z, 


dp(m)+ _ + . . . 

OZ z* 

(B 11) 

as z z c from above. It now follows from (B 6)-(B 8) that 


= o~ l M~ 2 — r ——{z ~ z c )-\ . 

* ma 0 c i z c 

(B 12) 

+ • • • , 

maoC! 9 z c 

(B 13) 

and 


ma 0 q z c 

(B 14) 


asz->2 c from above. It can be shown that the solution for is o(a~< M~ 2 ) as z -» 2 C 
from below, therefore the leading-order approximation to the velocity jump across the critical 
layer associated with the higher harmonics is given as 


Au 3 = u 3 (z+) - u 3 {z c ) = a 4 M 2 Re ^ 


771 = 2 


mao Ci 


6 a' m) e im ^ + 


(B 15) 
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Figure 6. Scaled instability-wave amplitude vs. long streamwise distance. 
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Figure 8(d). Vorticity contours in the (X, Y)-plane at x = 1.75. 
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